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The ground state properties and phase diagram of the bilayer square-lattice Heisenberg model 
are studied in a broad parameter space of intralayer exchange couplings, assuming an antifer- 
romagnetic coupling between constituent layers. In the classical limit, the model exhibits three 
phases: two of these are ordered phases specified by the ordering wave vectors (tt, tt; tt) and 
(0,0;7r), where the third component of each indicates antiferromagnetic orientation between 
layers, while another is a canted phase, stabilized by competing interactions. The effects of 
quantum fluctuations in the model with S — 1/2 have been explored by means of dimer mean- 
field theory, exact diagonalization of 2^2 x 2^2 x 2 clusters, and high-order perturbation 
expansions about the interlayer dimer limit. 
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§1. Introduction 

Quantum magnets in two dimensions have been a fascinating topic of numerous studies over the 
past decade. The interest has primarily focused on the competition between long-range magnetic 
order and novel disordered phases at T = as a result of enhanced quantum fluctuations. In 
this paper, we consider a bilayer square-lattice Heisenberg model with general nearest-neighbor 
exchange couplings in each layer and antiferromagnetic coupling between corresponding sites of 
each layer, described by the Hamiltonian 

H = Jl'^ Si_j-Sij + -^2 X! ^2,i-S2j + ^ Si,i-S2,i , (1.1) 

where the interlayer coupling is taken as the unit of energy. 

A particular limit of the model in which Ji = J2 = J > has attracted special attention 
for several reasons. Several cuprate superconductors contain Cu02 bilayers, so that model may 
be .eleva„t to thei. .a.„e«c p.ope„ieJ - and it is certainly a„ appropriate .odel fot the 
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magnetic properties of their parent compounds. Its phase diagram has been studied numerically 
by quantum Monte CarlJ) and perturbation expansions about J = 0,&'@^'& and analytically using 
linear spin-wave theory Ji> Schwinger boson mean-field theoryj^^ and an alternative bosonic mean- 
field calculation.^ This one-parameter subspace of the Hamiltonian (1.1) has two well known 
limits. For J ^ 1, the model describes a pair of two-dimensional Heisenberg antiferromagnets, 
weakly interacting with each other. In the isolation limit, each would be in the ordered state at 
T = and would possess Goldstone modes related to spontaneous breakdown of spin-rotational 
symmetry. Upon weak coupling, their sublattice magnetizations are phase-locked and half of the 
Goldstone modes acquire gaps. For J <^ 1, the ground state is very close to the "dimer" state, 
the product of spin singlets composed of adjacent interlayer pairs of spins. This state has a gap of 
order unity. Thus the system has a singlet ground state both for J ^ 1 and J <^ 1, and one can 
anticipate that there should be a critical value Jc of the intralayer coupling at which a continuous 
transition between the two phases is expected; furthermore, we can anticipate that this transition 
belongs to the universality class of the classical d = 3 Heisenberg model.& The various numerical 
studies all concur that Jc ~ 0.394. 

We were led to consider the more general Hamiltonian (1.1) on the following grounds. First of all, 
for J2 = the model is related to the "Kondo necklace" model introduced by Doniach0) as a simple 
model to describe the competition between the intra-atomic Kondo coupling (favoring local singlet 
states) and the inter-atomic RKKY interaction (favoring magnetic order and local high-spin states) 
in heavy fermion systems. The model for J2 = may also be viewed as the symmetric Hubbard- 
Kondo lattice model in the limit of strong on-site repulsion U between conduction electrons. B)© In 
the conventional symmetric Kondo lattice model (U = 0) on the square lattice with Kondo coupling 
Jk, it is known that a transition from spin-liquid to antiferromagnetic phase takes place when the 
transfer energy t is increased beyond tc ~ O.7Jft:.0^ Our calculation in §3.2 (the case 6 = 0) 



indicates that there is such a transition at tc ~ 0A2\/U Jk (for U ^ Jk)'- The antiferromagnetic 
order is suppressed by the strong on-site repulsion. This result is consistent with a finding by 
Shibata et a/.0^ for the one-dimensional symmetric Hubbard-Kondo lattice model, that the spin 
gap and inverse spin-spin correlation length are increasing functions of U. For J2 7^ 0, our results 
describe the effects of direct interactions between localized spins on the stability of the Kondo- 
singlet phase, which is neglected in the most analyses.0) As would be naively expected, the 
antiferromagnetic direct interactions are shown to favor the long-range ordered phase over Kondo 
singlets, while the opposite is true for (weak) ferromagnetic direct interactions. 

When Ji and J2 have opposite signs, the model is uniformly frustrated: The product of couplings 
around any four-spin plaquette that encompasses both layers is negative. Uniform frustration has 
been shown to lead to exotic magnetic behavior in many cases, both for theoretical mo dels© and 
real materials. Therefore it would be of great interest to see what are expected in the present 
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relatively simple model. We find no evidence for any exotic magnetic singlet phases, such as a 
spontaneously dimerized or plaquettized phase. 

We have carried out several types of analyses in order to identify the phases and phase transitions 
exhibited by the Hamiltonian (1.1). The classical ground state phase diagram is first considered 
(in §2). Then quantum effects for S = 1/2 spins are included, and it is clearly shown that a 
spin-disordered phase (the aforementioned dimer phase) appears in the vicinity of Ji = J2 = 0. 
We examine the instabilities of this phase by means of mean-field theory, small-system exact diag- 
onalization (in §3.1), and higher-order perturbation expansions for susceptibilities and excitation 
spectra (in §3.2). The last technique also allows us to determine the spin- wave velocity in the 
magnetically ordered states in the vicinity of the phase boundary. The last section (§4) is devoted 
to summary and discussion. In particular, Sachdev and Senthil0) recently considered a quantum 
rotor model which is closely related to the model (1.1), and in §4 we will compare their results with 
ours. 

§2. The classical limit 

For our purpose in the following, it is instructive to consider the ground state phase diagram in 
the classical limit, where the operators S in the Hamiltonian (1.1) are replaced by classical vectors. 
By noting that the interchange Ji J2 simply corresponds to relabeling the layers, we only need 
to consider the half space Ji > J2 of the J1-J2 plane. 

One can easily determine the ground states in the part of parameter space where the two intralayer 
couplings have the same sign. When Ji and J2 are both positive the ground state is composed of 
a pair of Neel ordered layers with an ordering wave vector (7r,7r;7r), where the third component 
indicates the antiferromagnetic orientation between layers. When Ji and J2 are both negative 
the ground state is composed of two ferromagnetically (but oppositely) ordered layers with an 
ordering wave vector (0, 0;7r). In each case, both intra- and interlayer energies of the system can 
be minimized simultaneously. 

The (7r,7r;7r) phase and (0,0; vr) phase both extend into the quadrant Ji > > J2, bordered on 
a phase which will be referred to as the weakly ferromagnetic phase (WF). This phase is composed 
of a canted Neel-ordered and a canted ferromagnetic layers, as is shown in Fig. 1. The direction of 
the staggered moment in the antiferromagnetically coupled layer is perpendicular to the uniform 
moment on the ferromagnetic layer. The canting angles ip and (j) in the layers satisfy 




(2.1) 



and 




(2.2) 



with defined by 




(2.3) 
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When Ji and J2 are large, the canting angles behave like ^ l/(8Ji) and (/) 1/(8J2). 

The boundary between the (0,0; tt) phase and the WF phase is the line 1/Ji + I/J2 = 8, while 
the boundary between the (tt, 7r;7r) phase and the WF phase is the line 1/Ji + I/J2 = —8. In both 
cases the phase transition is continuous. The classical phase diagram is plotted in Fig. 2. 



Fig. 1. Canting angles ^p, (p in the antiferromagnetic (AF) and ferromagnetic (F) layers respec- 
tively. 



Fig. 2. Phase diagram in the classical spin limit {S — > 00). 
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§3. Phase diagram for the quantum model 

3.1 Simple considerations 

There is at least one region in parameter space where the classical phase diagram, Fig. 2, should 
be modified for the quantum spin model with S = 1/2. When Ji and J2 are both small, the ground 
state is approximately a product of interlayer-dimer singlets and does not exhibit long-range order 
in its spin-spin correlation function. Since there is an energy gap of 1 (the interlayer coupling) for 
Ji = J2 = 0, we can expect that the "dimer" phase is stable in a region enclosing the origin with 
diameter of order unity. Let us try to make this statement more precise, by means of relatively 
simple calculations. 

An estimate of the boundaries between the dimer phase and long-range ordered phases [{it, tt; tt), 
(0, 0;7r), WF] can be obtained from dimer mean-field theory. One considers the following two-site 
Hamiltonian 

Si • S2 - hiSf - h2Sl - g^Sf - g^S^ (3.1) 
and imposes the self-consistency conditions 

hi = AJi{Sl) h2 = U2{SI) gi=g2 = (3.2) 

for the (tt, tt; tt) phase, or 

hi = -4Ji{Sf) h2 = -4J2{SI) 51 = 52 = (3.3) 

for the (0, 0; tt) phase, or 

hi=AMSf) /i2 = 4J2(S|) gi = -4Ji{Sf) 52 = -4J2(-Sf) (3.4) 

for the WF phase. In the above self-consistency conditions, the angular brackets denote ground- 
state averages. 

As usual in mean-field theory, one finds that /ii = /12 = 51 = 52 = is the only solution, 
and hence the dimer phase is stable, for a small inter dimer couplings. It is straightforwardly 
shown that the dimer phase is unstable against the (7r,7r;7r) phase for Ji + J2 > 1/2, and the 
(0,0; tt) phase for Ji + J2 < —1/2; in both cases the transitions arc continuous. Much more 
effort is needed to determine the boundary of the WF phase, and we have not been able to derive 
any analytical results. Instead, we have carried out numerical calculation, solving the mean-field 
equations iteratively (starting with a set of prescribed magnetizations to calculate the fields and 
solving the two-site Hamiltonian for the magnetizations, repeating these process until convergent 
results are obtained). The resulting phase boundaries are shown in Fig. 3. An unphysical feature 
of our results is that the ground-state energy exhibits a discontinuous decrease on entering the WF 
phase. This might be due to the WF solutions which are not stable under iteration, so that we 
have underestimated the domain of stability of the WF phase. 
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We have also considered another type of mean-field calculation, which was applied to the special 
case Ji = J2 by Chubukov and Morr.i^ Their calculation can be readily generalized. Let us omit 
all details, which are clearly described in their paper, and just describe the results. By rewriting 
the Hamiltonian in terms of three types of bosonic operators, corresponding to the three triplet 
excited states of an 5* = 1/2 dimer, and keeping only the quadratic terms, one ends up with a 
three-fold degenerate excitation spectrum of the form 

e(q) = ^1 + (Ji + J2)icosqx + cosqy) . (3.5) 

(The in-plane lattice spacing is taken as the unit of length.) At small intralayer couplings the 
spectrum has a gap, characteristic of the dimer phase. When | Ji + J2I is increased beyond 1/2, 
some of the excitation energies pass through zero onto the imaginary axis, signaling an instability 
of the dimer phase. This analysis agrees completely with the dimer mean-field calculations for the 
transitions to the (0, 0; vr) and (vr, vr; it) phases, as described above. It also yields the result that 
the spin-wave velocity on both of those phase boundaries is given simply by c = 1/2. However, this 
analysis does not appear capable of describing the transition to the WF phase. 



Fig. 3. Phase diagram for 5=1/2 spins from dimer mean-field theory. 

Finally, we have carried out finite-size exact diagonalization studies on 2 x 2 x 2 and 2\/2 x 2-v/2 x 2 
systems with periodic boundary conditions.0^ In Fig. 4 the low-lying energy levels are plotted as 
functions of J = \ Ji\ = | J2I for the three distinct choices of the signs of couplings. One finds that 
when the intralayer couplings have the same sign the ground state is always a singlet and when the 
couplings have opposite signs there are level crossings to high-spin states as the magnitude of the 
couplings increases. The latter is associated with the transition from the dimer phase to the WF 
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phase; in the WF phase the ground state has a net moment and therefore must have nonzero spin. 
We will have more to say about the finite size calculations in the following subsection. 



Fig. 4. Low- lying energy levels of 2 x 2 x 2 (triangles) and 2\/2x 2\/2 x 2 (dots) clusters as functions 
of J = I Ji| = I J2I : (a) Ji > 0, J2 > 0, (b) Ji < 0, J2 < 0, (c) Ji > 0, J2 < 0. 

3.2 Series expansions and extrapolations 

In order to obtain more quantitative estimates of the domain of stability of dimer phase and the 
spin-wave velocity in the magnetically ordered phases adjacent to those phase boundaries, we have 
carried out perturbation expansions in powers of the intralayer couplings about the pure interlayer 
dimer Hamiltonian. 

For convenience, the intralayer couplings are parametrized by 

Ji = A cos 9 J2 = A sin ^ . (3.6) 
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Expansions in A can then be carried out at fixed 6, rather than carrying out a double expansion 
in powers of Ji and J2 simultaneously. Note that only expansion for 9 in the range — 7r/4 to 7r/4 
is needed, since negative values of A correspond to 6 in the range 37r/4 to 57r/4, which then maps 
to —it/ 4 to — 37r/4 under the interchange Ji <-> J2. In this way full 27r coverage is obtained by 
carrying out calculations only for a quarter of the circle, which were done at intervals of 7r/16. 
Series expansions for several susceptibilities xi^) at T = and the Fourier transform of the triplet 
elementary excitation spectrum e(q) were performed up to the order A^ using connected cluster 
methods. 

The values of A corresponding to continuous transitions out of the dimer phase can be accurately 
estimated by constructing differential approximants0^ to the susceptibility and the gap series asso- 
ciated with the ordering wave vector. The leading critical behavior of these quantities is described 
by X ~ i^c — A)""^ and e ~ (Ac — A)'^. The dimer-(7r, tt; tt) and dimer-(0, 0; vr) transitions should lie 
in the universality class of the d = 3 classical Heisenberg model, and we can expect 7 ~ 1.4 and 
along those critical lines. If the dimer- WF transition is continuous, the naive hypothesis 
is that it would lie in the d = 3 classical 5*0 (3) class, as has shown to be the case for canted magnets 
by Kawamura, i) for which 7 ~ 1.1 and u ~ 0.53. On the other hand, if the dimer- WF transition 
is of the first-order, one would expect the linearly vanishing gap, since the transition would be a 
result of a level crossing to a state of different symmetry. The value of A at which the gap vanishes 
would be regarded as an upper bound on the actual phase boundary, since the excitations might 
have attractive interactions (and as we will see below, this is the case in the relevant parameter 
regime) . 

In order to organize the following discussion, it is useful to divide the parameter space into three 
regions: region I, with < 9 < ir/4:; region II, with — 37r/4 < 9 < —it/2, and region III, with 
-tt/2 <9 <0 . 

In region I, the dimer phase is adjacent to the (vr, tt; tt) phase, and a continuous transition between 
these phases is expected as in the conventional bilayer Heisenberg antiferromagnet {6 = tt/4). 
Shown in Fig. 5 is the phase diagram based on biased differential approximants to the expansion 
series x(7'")7r;7r) and e(7r,7r). 

Note that without biasing, the critical exponents derived from approximants show good agree- 
ment with the anticipated values, but the biasing allows for more precise estimates of critical 
couplings than would be possible otherwise. The uncertainties in the critical value for A estimates 
are approximately 0.002. 

In region II, there is a continuous transition from the dimer phase to the (0, 0; vr) phase, and 
again the critical exponents 7 and u agree well with the anticipated values. Critical points derived 
from biased differential approximants are shown in Fig. 5. 

It is to be noted that the inequality | Ji + J2I > 1/2 along these two lines of continuous transitions 
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is held, indicating that the mean-field theory underestimates the domain of stability of the dimer 
phase. This is entirely plausible, since mean-field theory neglects the effects of quantum fluctua- 
tions, which stabilize the dimer phase. Furthermore we note that these two lines are not mirror 
images of one another. The dimer-(7r, vr; vr) line is slightly but noticeably curved. The curvature 
suggests that quantum fluctuations suppress Neel order more strongly when the coupling between 
dimers is evenly distributed between the layers than when the coupling is concentrated in one layer. 
On the other hand, the critical dimer-(0, 0; vr) line is almost as straight as is the case in dimer mean- 
field theory. We have not an clear explanation for it, but it is worth noting another point which 
shows that the (vr, vr; vr) and (0, 0; vr) phases are affected differently by the quantum fluctuations: 
Within linear spin-wave theory, the spin-wave velocity is given by [(Ji + J2)/4 + 2Ji J2]i/2 in the 



represent precise estimates of continuous transitions; the dashed lines are more uncertain estimates 
of transitions, whose characters are discussed in the text. The dot-dashed lines are based on finite- 
size calculations. 

In region III, the situation is more complicated than in the preceding two. Based on the classical 
phase diagram and dimer mean-field theory, one expects that the dimer-(vr, vr; vr) and dimer-(0, 0; vr) 
critical lines would extend somewhat into this region, and this is in fact observed. On increasing 




Fig. 5. 



Phase diagram for S — \/2 spins based primarily on series expansions. Solid lines 
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A, a dimer-(7r, tt; tt) transition takes place at ^ = — 7r/16, and likewise a dimer-(0, 0; tt) transition 
occurs for = — 77r/16: the relevant gaps vanish and susceptibilities diverge at the same value of A 
with the expected exponents. 

However, for most values of 9 one would anticipate encountering dimer-WF phase transitions. 
In fact there is an evidence for such transitions from unbiased approximants to the series, but 
the detailed features are rather complicated, in that neither the first-order nor continuous phase 
transition scenario receives complete support. Along the dashed line in Fig. 5, one finds that e(7r, tt) 
vanishes with exponents close to 1/2. The staggered susceptibility diverges with exponents in the 
range 0.8-1.0. So far, this is consistent with the continuous transition scenario. However, the 
critical values of A for the susceptibility are typically 20% larger than those for the gap, and the 
"error bars" (deduced from the consistency of the unbiased approximants) of two sets of critical 
points do not overlap with each other. Since the exponents of susceptibility are generally smaller 
than theoretically anticipated, one might hope that a biased analysis of the susceptibility would 
bring the two sets of estimated critical points into agreement. Unfortunately, this is not the case 
and biasing actually drives them farther apart. Consequently, we cannot rule out the possibility 
that the dimer-WF transition involves a vanishing gap but finite susceptibility, that is, a first-order 
transition. 



Fig. 6. A map of the total spin of ground state of 2\/2x2\/2x2 clusters in the J1-J2 plane. 

As for the phase boundaries between the WF and magnetically ordered phases, the series expan- 
sions provide no information. The mean-field and classical calculations suggest that they lie nearly 
parallel to the Ji and J2 axes. In order to obtain a crude estimates of those phase boundaries, we 
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have constructed a map of the total spin of ground state of 2\/2 x x 2 clusters in the J1-J2 
plane, determined by exact diagonalization, see Fig. 6. The region where the ground state has 
nonzero spin is presumed to lie in the WF phase. The dimer-WF boundary estimated in this way 
is consistent with that derived from the series analysis. Therefore we expect that the boundaries 
between the WF and magnetically ordered phases (shown by the dot-dashed lines in Fig. 5) are 
reasonably accurate. 

At this point let us consider the properties of the excited states of the 2\/2 x 2\/2 x 2 clusters 
in more detail. We find that in the vicinity of the singlet-singlet (dimer-to-antiferromagnetically 
ordered) phase transitions the lowest excitation energy in the 5 = 2 sector is twice greater than 
that in the 5 = 1 sector, while in the vicinity of the dimer-WF phase boundary the reverse is true: 
see Fig. 7. 



Fig. 7. Excitation gaps (A = £;nim(5') - E^i^(S = 0)) scaled by S or 5(5 + 1) in 2^2x2^2x2 
clusters as functions of A: (a) Q = tt/4, (b) 6 = — 7r/4. 

Thus the triplet excitations are repulsive in the former cases and attractive in the latter, as antici- 
pated. We have also examined how the lowest excitation energies in each spin sector vary with S 
in each of the singlet ground states. In the spin-disordered dimer phase those energies are roughly 
proportional to S, while in the ordered phases they are nearly proportional to S{S + 1). Both 
results are readily understood. In the dimer phase the lowest spin-S state is essentially composed 
of S triplet excitations. In the magnetically ordered phases the situation is rather different. In 
the bulk, arbitrarily low-energy excitations of any size of spins are possible, but for large but finite 
clusters the excitations amount to finite angular momentum states of the staggered moment (for 
which the corresponding moment of inertia is the uniform susceptibility) which accounts for the 
S{S + 1) scaling.!) The crossover between the two different behaviors is impossible to pin down 
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the critical points precisely in a system as small as the one we have considered, but the critical 
points obtained from the series analysis are clearly in the range where the crossover is taking place, 
see Fig. 7 (a). 

Finally let us present further results obtained from the series expansions regarding properties 
along the boundary of the dimer phase. The spin-wave velocity on the dimer-(7r, vr; vr) and dimer- 
(0,0; vr) critical lines was determined using the method described in Ref. ^). In Fig. 8, we show 
how the quantum renormalization factor of the critical spin-wave velocity defined by 

Zc(^) = c/cLswT (3.7) 

depends on 0, where clswt is the spin- wave velocity within linear spin- wave theory whose value 
was given earlier. To a good approximation we find that Zc{0) is independent of 6 along the dimer- 
(0, 0; vr) boundary. This seems to be consistent with the fact that the phase boundary itself can be 
well described by Ji + J2 = const. 



Fig. 8. Quantum renormalization factor of the critical spin- wave velocity Zs{9) along continuous 
transition boundaries as a function of 9. 

One can also extrapolate the excitation spectra for values of q far from the characteristic mo- 
menta (0,0) or (vr,7r), where one expects direct Fade approximants to provide accurate estimates. 
Estimates of the dispersion relations along selected lines in the Brillouin zone are presented for 
several 6 values in Fig. 9. 
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Fig. 9. Critical spin-wave excitation spectra: solid, dotted and dashed lines are for = 0, 
6 = —5n/16 and 6 = — 7r/2, respectively. 

§4. Summary and Discussion 

We have investigated the ground state properties and phase diagram of bilayer Heisenberg model 
on the square-lattice with general nearest-neighbor intralayer couplings. The ground state phase 
diagram of the model has been first examined in the classical limit in order to get an insight to its 
asymptotic behaviors. Then the effects of quantum fluctuations are studied on the basis of several 
methods; dimer mean-field theory, small cluster diagonalization, and high-order perturbation the- 
ory about the interlayer dimer limit. We have found four phases for the model with S = 1/2: one 
"spin liquid" (dimer) phase which is the analytic continuation of the interlayer dimer Hamiltonian, 
two antiferromagnetically ordered [(vr, vr; vr) and (0, 0; tt)] phases, and one phase with weak mag- 
netization stabilized by competing interactions between layers (WF). The transitions between the 
dimer phase and the antiferromagnetically ordered phases are continuous (and belong to the clas- 
sical d = 3 Heisenberg model universality class), while the character of the dimer- WF transitions 
could not be reliably determined. Concerning the dimer-(7r, vr; vr) and dimer-(0, 0; vr) transitions, it 
appears that quantum fluctuations play a more important role in the former, and lead to a rather 
intriguing result that Neel order is suppressed more strongly when the intralayer couplings are 
evenly distributed between the layers than when they are concentrated in one layer. 

The case {9 = 0) is especially interesting. In this case, the spins are antiferromagnetically 
coupled in one of layers and not coupled at all in the other. This model corresponds to symmetric 
Hubbard-Kondo square-lattice model in the large-C/ limit. Provided t ^ oo and U co with 
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At'^/U = Js = constant (the intralayer coupling) and the on-site Kondo couphng (the inter layer 
coupling) is Jk, we have found that there is a transition to a magnetically ordered state when Jg 
is increased up to approximately 0.71Jk- This leads to the result for the critical t mentioned in 
the Introduction. 

The one-dimensional version of the Hamiltonian (1.1), with general intrachain couplings, has 
been attracting interest recently in connection with the spin-ladder materials® In particular, 
Tsukano and Takahashi@) studied the one-dimensional model with J2 = —Ji {0 = — vr/4) and found 
a continuous transition between the dimer phase and a magnetically ordered phase, supporting 
continuous transition scenario of the dimer- WF transition. 

Finally, Sachdev and Senthil0) have studied a rather general quantum rotor model for which, 
in the large-S limit, there is a direct mapping to the model (1.1). It is interesting to see which 
features of their model persist for S = 1/2, and which do not. The mean-field analysis of the 
rotor model finds a "gapped quantum paramagnet," equivalent to the dimer phase in the spin 
model, a Neel phase, a "quantized ferromagnet" phase, and a canted phase, equivalent to the WF 
phase in the spin model. They suggest that, quite generally, continuous phase transitions into the 
canted phase, where there is nonzero magnetization, are only possible from phases with gapless 
excitations. Consequently, continuous phase transitions between the dimer phase and the WF 
phase should not be possible. (In the section of the phase diagram for the rotor model which is 
presented by Sachdev and Senthil, not even first order transitions between the gapped quantum 
paramagnet and canted phases are exhibited, but there is no reason to believe that holds for all 
choices of parameters.©) Our numerical studies were not able to establish the nature of the dimer- 
WF transition, and so our results are not in direct disagreement with the general principle espoused 
by Sachdev and Senthil even though they do not provide any support for it, either. However, there 
is one significant difference between the phase diagrams of the rotor model and of the spin model. 
For S = 1/2, the Hamiltonian (1.1) does not exhibit an analog of the quantized ferromagnet phase. 
Ferromagnetically coupled bilayers of S = 1/2 spins could exhibit such a phase, since in the absence 
of intralayer coupling the interlayer dimers would have angular momentum 1. In contrast to the 
rotor model, however, it does not appear possible to go to such a phase directly from the dimer 
phase by varying the intralayer couplings. 
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